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, Abstract 

Motivated by the immense success of random walk and Markov chain methods in the design of 
classical algorithms, we consider quantum walks on graphs. We analyse in detail the behaviour of 
unbiased quantum walk on the line, with the example of a typical walk, the "Hadamard walk" . We 
show that after t time steps, the probability distribution on the line induced by the Hadamard walk is 
almost uniformly distributed over the interval [—t/s/2, t/y/2]. This implies that the same walk defined 
on the circle mixes in linear time. This is in direct contrast with the quadratic mixing time for the 
corresponding classical walk. We conclude by indicating how our techniques may be applied to more 
general graphs. 

Oh 



1 Introduction 

Random walks on graphs have found many applications in computer science, including randomised al- 
gorithms for 2-Satisfiability, Graph Connectivity and probability amplification (see, e.g., |14j|). Recently, 
Schoning |r| discovered a random walk based algorithm similar to that of [17] that gives an elegant (and 



the most efficient known) solution to 3- Satisfiability. In general, Markov chain simulation has emerged as 
a powerful algorithmic tool and has had a profound impact on random sampling and approximate count- 



ing [10]. Notable among its numerous applications are estimating the volume of convex bodies ||[] and 
approximating the permanent ||. A few months ago, Jerrum, Sinclair and Vigoda [|ll|] used this approach 
to solve the long standing open problem of approximating the permanent in the general case. 

In the spirit of developing similar techniques for quantum algorithms, we consider quantum walk on graphs. 
To date, few general techniques are known for developing and analysing quantum algorithms: Fourier 



sampling, which is typified by the seminal work of Simon [21] and Shor pCfl , and amplitude amplification, 
which originated in the seminal work of Grover |8|. Barring applications of these techniques, the search for 
new quantum algorithms has primarily been ad hoc. We believe that studying quantum walk on graphs is 
a step towards providing a systematic way of speeding up classical algorithms based on random walk. 

* Supported by a joint DIMACS- AT&T Post-Doctoral Fellowship and NSF grant EIA 00-80234. 

' Supported by a C.E. Procter Fellowship. 

1 See [|L2j for the latest progress on this problem. 
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Figure 1: The dynamics of the Hadamard walk. In (a) we begin at time t with a particle in chirality state 
"right" or "left". The result of the Hadamard transformation is shown in (b), the particle is now in an 
equal superposition of "left" and "right" chirality states (with the amplitudes indicated) and then moves 
accordingly (c) to generate the state at time t + 1. 



We begin by considering the quantum walk in one dimension. Recall that a classical particle doing a 
random walk on the integer lattice chooses, at every time step, a random direction ("left" or "right") to 
move in, and moves to the site adjacent to it in that direction. In direct analogy, one may naively try to 
define quantum walk on the line as follows: at every time step, the particle moves, in superposition, both 
left and right with equal amplitudes (perhaps with a relative phase difference). However, such a walk is 
physically impossible, since the global process is non-unitary. It is also easy to verify that the only possible 
homogeneous (i.e., translationally invariant) unitary processes on the line involving transitions between 
adjacent lattice sites are the left and right shift operators (up to an overall phase) |{13[| . This corresponds 
to the rather uninteresting motion in a single direction. As explained below (and also shown by fl3|| ), it is 
still possible to construct a unitary walk if the particle has an extra degree of freedom that assists in its 
motion. 

Consider a quantum particle that moves freely on the integer points on the line, and has an additional 
degree of freedom, its chirality, which takes values "left" and "right" . A walk on the line by such a particle 
may be described as follows. At every time step, its chirality undergoes a rotation (a unitary transformation 
in general) and then the particle moves according to its final chirality state. Figure [l] depicts this two-stage 
move in a quantum walk where the chirality undergoes a Hadamard transformation. We call this particular 
walk the "Hadamard walk". 

In this paper, we analyse in detail the dynamics of the Hadamard walk on the line. We derive the asymptotic 
form of the probability induced on the line by observing the position of a particle doing the walk and show 
that after t time steps, the distribution is almost uniformly distributed over the interval [-t/y/2, t/y/2}. 
This implies that the analogously defined walk on the or the circle mixes in linear time. This is in 
immediate contrast with the classical random walk, which mixes in quadratic timeF] 

2 In this paper, by the 5-mixing time of a walk we mean the first time (independent of the initial state) at which the 
distribution induced by the walk is 5-close to uniform in total variation distance. The parameter 5 is understood to be a 
constant less than 1 if it is not explicitly specified. 
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Next, we turn to general quantum walk on the line. Such a walk may involve arbitrary transitions between 
adjacent lattice points on the line (subject to the global unitarity constraint). It turns out that any 
homogeneous, local unitary process on the line governing the motion of a particle may be recast as a 



unitary transformation of the chirality followed by a left/right shift [16|. Moreover, if we are interested in 
unbiased walk[], we may assume that the particle moves one step to the left if it ends up with chirality 
"left" , and one step to the right if it ends up with chirality "right" . Thus, there is no loss in generality in 
studying only such walk. 

It can be shown fl6fl that if we are interested only in the probability distribution induced on the line by such 
a walk, we may restrict ourselves to the study of a family of walks specified by a single parameter 9 £ [0, tt}. 
We show that every walk in this family, except that for 6 = 0, it, shares the gross features of the Hadamard 
walk. In particular, all except the two singular walks mix in linear time. 

One method for analysing quantum processes such as the walk we study is the path integral approach. 
This method involves explicitly computing the amplitude of being in a certain state as the sum, over all 
possible paths leading to that state, of the amplitude of taking that path. The problem of evaluating, or 
even approximating, path integrals to determine the time development of a wave function is well known to 
be notoriously hard. We circumvent this barrier by taking the Schodinger approach, which takes advantage 
of the time and space homogeneity of the quantum walk. The crucial observation is that because of its 
translational invariance, the walk has a simple description in Fourier space. The Fourier transform of the 
wave function is thus easily analysed, and transformed back to the spatial domain. It is noteworthy that 
this technique is standard in the analysis of classical random walk || . 

A key advantage of the Schodinger approach is that the resulting description of the wave function, in terms 
of integrals of the Fourier type, is amenable to analysis in standard ways. There is a well-developed theory 
of the asymptotic expansion of integrals Q which allows us to determine the leading behaviour of the wave 
function in the large time limit. (This is exactly the case of interest in the design of asymptotically efficient 
algorithms.) We are thus able to derive the asymptotic functional form of the corresponding probability 
distribution. 

We also use the theory of asymptotic expansion of integrals in a non-standard way to deduce some properties 
of the probability distribution induced by the walk. For example, we use this method to approximate the 
probability masses of each point on the line by an integral of the Fourier type. This allows us to calculate 
the net probability of an interval, and the different moments of the distribution very precisely. 

Finally, we note that the asymptotic quantities we calculate match very well with simulation results even 
for small times, indicating that our results may apply to small quantum systems as well. 



Related work 



Quantum dynamics similar to ours has previously been studied by Meyer [13] in the context of quantum 
cellular automata. He singles out the one parameter family of unitary processes mentioned above and 
derives a closed form solution for the corresponding wave functions. The solution is obtained by following 
the path integral approach and allows him to express it in terms of the Gauss hypergeometric function. His 
motivation in studying such dynamics is however very different from ours, and he does not quantitatively 
analyse the wave functions further. 



3 Appendix [d| has a formal definition of this notion. Essentially, it means that the walk makes no distinction between the 
two halves of the line given by its initial position. It is not hard to see that in any other kind of (biased) unitary walk, the 
particle moves in a single direction with probability 1 at every time step. 
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In work more closely related to ours in its motivation, Nayak, Schulman and Vazirani [15] have studied 



random sampling on a quantum computer. They gave an efficient method for generating superpositions for 
which the corresponding probability distributions satisfy certain well-behavedness properties. Using this 
for a distribution related by a Fourier transform, they gave an efficient quantum algorithm for sampling 



from the Gibbs distribution for the Ising model. (Subsequently, Randall and Wilson [18] discovered a 
classical polynomial time algorithm for this task via a different route.) 



Watrous (2^] has considered unitary processes based on quantum walk on regular graphs in the context 
of Undirected Graph Connectivity and logarithmic-space computation. He shows that it is possible to 
construct, in logarithmic space with limited measurement, and with high probability, a good approximation 
to the superposition over all vertices connected to a specified vertex in the graph, using a local unitary 
transform subroutine. 

We have recently learnt that Aharonov, Ambainis, Kempe and Vazirani [jlj] have independently studied 
quantum walks on graphs and their mixing behaviour. They observe that the probability distribution 
induced by a unitary process on a finite dimensional Hilbert space is essentially periodic and hence does 
not converge to a stationary distribution. They argue that the notion of mixing in this case is therefore 
more appropriately captured by the closeness to the uniform distribution, of the average of the probability 
distribution over all time steps. They show that the quantum walk on the n-circle mixes in time 0(n log n) 
in this sense. They also show a lower bound of 1/cM? for this mixing time for general graphs, where d is 
the maximum degree of a vertex in the graph and <3? is the conductance. 

We have also learnt that Ambainis and Watrous have analysed the Hadamard walk on the unbounded 
line, and show that it has linear mixing characteristics. They follow the approach of Meyer, and further 
analyse the path integrals so obtained to arrive at the result. 



Organisation of the paper 

The rest of the paper is organised as follows. We formally define and analyse quantum walk on the line in 
Section ||. Some background on Fourier transforms required for this is summarised in Appendix [A]. Details 
of some calculations required for the section are provided in Appendix |B|. We then consider the behaviour 
of the walk after large times and derive its asymptotic properties in Section |3[ This section is heavily based 
on the Method of Stationary Phase described in Appendix [A[ The approximations made in the section are 
justified in Appendix ^|. For lack of space, we are not able to provide the full details of the analysis of the 
general quantum walk on the line. The technique involved in the analysis is however the same as that for 
the Hadamard walk. Appendix ^ summarises the main conclusions about the general case. 



2 One dimensional quantum walk 
2.1 Formal description of the walk 

A particle doing a (classical) symmetric random walk on the line may be described as follows. It starts, 
say, at the origin, and at every time step, tosses a fair coin. Each of the two possible outcomes of the 
toss is associated with a distinct direction, "left" or "right" . The particle moves one step in the direction 
resulting from the toss. 

A quantum generalisation of this process involves a quantum particle on the line, but with an additional 
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degree of freedom which we call the "chirality" . The chirality takes values "left" and "right" , and directs 
the motion of the particle. At any given time the particle may be in a superposition of "left" and "right" 
chirality states and is therefore described by a two-component wave function. The dynamics of the walk 
that we consider is given by the following rules. At each time step, the chirality undergoes a rotation (a 
unitary transformation, in general), and the particle moves according to its final chirality state. Therefore, 
if the particle ends up with chirality "left", it moves one step to the left, and if it ends up with chirality 
"right", it moves one step to the right. 

For concreteness, we begin by focusing on a quantum walk in which the unitary transformation acting on 
the chirality state at each time step, is chosen to be the Hadamard transformation (the Fourier transform 
over Z2): 

|L) - -±=(|L) + |R» 
|R) - ^(|L)-|R». 

Here L and R refer to the "right' and "left" chirality states. The resulting quantum walk, which we will 
refer to as the "Hadamard walk" is depicted in Figure ||. 

To study the properties of the walk defined above, we consider the wave function describing the position of 

ih(n,t) 
^ ipR{n,t) 

of amplitudes of the particle being at point n at time t, with the chirality being left (upper component) or 
right (lower component). The dynamics for is then given by the following transformation (cf. Figure |l|): 



the particle and analyse how it evolves with time. Let ^(n, t) 



be the two component vector 



*(n,i + l) 



*(n-l,t) + 





1 -1 

M + V(n - l,t) +M_*(n + l,t) 



1 1 




*(n + l,i) 



for matrices M + ,M_ defined appropriately. Note that this transformation is unitary on the basis states 
given by Z x since it is the composition of a unitary operation and a reversible move to the left or the 
right. Moreover, since the particle starts at the origin with chirality state "left" (say), we have the initial 

I ), and¥(n,0) = ( ° Q 



conditions, ^(0,0) 



if n ^ 0. 



With the above formulation, the analysis of the Hadamard walk reduces to solving a two dimensional linear 
recurrence (or a difference equation). We show how this recurrence may be analysed in the next section. 



2.2 Fourier analysis of the Hadamard walk 

As mentioned in Section [l], quantum walk of the kind defined above has, due to translational invariance, 
a very simple description in the Fourier domain. We therefore cast the problem of time evolution in this 
basis, where it can be easily solved, and at the end revert back to the real space description by inverting 
the Fourier transformation. The details are described below. 

The spatial Fourier transform ty(k,t) (for k € [— 7r, 7t]) of the wave function ty(n,t) over Z is given by (cf. 
Appendix [A]) 

¥(M) = £*(rM)e ifc ". 

n 
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In particular, we have ^(k,0) = y J J for all k, for a particle starting at the origin with initially "left" 
chirality. 

From the dynamics of ^, we may deduce the following about \E': 

*(M + i) = ^(M + *(n-i,t)+M_tf(n + i,t)) eikn 

n 

= e ik M+ *(" ~ M) e ifc ( n -^ + e~ ifc M_ ^ *(n + 1, t) e ifc ( n+1 ) 

n n 

= (e ifc M + + e~ ifc M_) *(jfe, *). 



Thus, we have, 
where 



*(M + i) = M**(M) 



M fc = e ifc M+ + e" ifc M„ (2) 
1 — 0) 



v/2 



e -ifc e -ifc 



Note that = A^?7 T , where is the diagonal matrix with entries e _lfc ,e lfc and U T is the transpose of 
the unitary matrix U that acts on the chirality state of the particle. Hence M/% is a unitary matrix. [This 
general presentation anticipates the walk with general U.] 

The recurrence in Fourier space thus takes the simple form *f>(k,t + 1) = *$>(k,t), leading to fy(k,t) = 
M|\I/(A;, 0). We may calculate Mt (and thus ^f(k,t)) by diagonalising the matrix M/%, which is readily 
done since it is a 2 x 2 unitary matrix. If has eigenvectors (|3>jt),|$!)) an( ^ corresponding eigenvalues 

M k = \l\<S>l)(<S>l\ + Al|$i)($l|, 
and then immediately we obtain the time evolution matrix as: 

The eigenvalues of are \\ = e~ luJh and A 2 , = e x ^ +ulk \ where ujk is defined as the angle in [— f , f] such 
that sin(cjfc) = The corresponding eigenvectors are displayed in Appendix [B|. For definiteness, here 

we consider the time evolution of a particle that begins at the origin in the "left" chirality state. In the 
Fourier basis this initial state is represented by Sfr(k,0) = (1,0) T for all k. The wave function at time t is 
then given by: 

r , . 1 COS k , • , ( — !)*/. COS/C , , s 

i/j L (k,t) = -(! + —= )e~"" fc * + v ) e"^* (4) 

2 v 1 + cos 2 2 v 1 + cos 2 k 

= - (e-^-(-l)V^) (5) 

2y 1 + cos z k 

We now invert the Fourier transformation, to return to the basis in real space (cf. Appendix [A|). The wave 
functions in real space can be written in the form: 
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Probability Distribution for Hadamard Walk:Asymptotics vs Simulations 



Asymptotic 
Simulation 




Figure 2: A comparison of two probability distributions, one obtained from a computer simulation of the 
Hadamard walk, and the other from from an asymptotic analysis of the walk. The number of steps in the 
walk was taken to be 100. Only the probability at the even points is plotted, since the odd points have 
probability zero. 



Mn , t) = LUdTl I" * (1 + , (6 , 

I J--* 2ir a/1 + cos 2 k 

Mn . t) = i±^r*_*^.~> (7) 

1 J-tt 2ir a/1 + cos 2 k 

Notice that the amplitudes vanish for even n (respectively, odd n) at odd t (even t), as we would expect 
from the definition of the walk. 

Thus, we have obtained a closed form solution for the time evolved wave function of the Hadamard walk. 
In view of possible applications to developing Quantum algorithms, we are naturally led to considering the 
behaviour of the wave functions at large times (t > 1). Happily, the problem is considerably simplified in 
this asymptotic limit which allows us to accurately derive several useful results. In the next section, we 
give details of this asymptotic analysis. 



3 Asymptotic properties of the wave function in the large time limit 

In the previous section, we obtained an exact solution for the time evolution of the Hadamard walk. In 
what follows, we will use extensively the Method of Stationary Phase (see Appendix |A|) to extract the 
asymptotic properties of the resulting wave function. 
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3.1 The asymptotic probability distribution 

The asymptotic analysis for tpi, and V'R is essentially the same. They can both be written as a sum of 
integrals of the type I(a,t) described in Appendix [B|. Here, a should be understood as n/t. Appendix |B| 
gives the details of the analysis for this generic integral; we summarise the conclusions below. 

The asymptotics for the wave function are simple to describe. The wave function is essentially uniformly 
spread over the interval between =pl/\/2 (and so there its gross behaviour is like 1/y/i), and it dies out very 
fast (faster than any inverse polynomial in t) outside the interval. At the "frontiers" at =pl/\/2, however, 
there are two peaks of width Oft 1 / 3 ) where the wave function goes as i -1 / 3 . The net probability below 
these two peaks and beyond =Fl/\/2 thus goes as i -1 / 3 , and does not affect the properties of the wave 
function we are interested in. We will therefore restrict our attention to the interval a E (— 1/\/2, \j\pl). 
The exact expressions for ^ = (V'L)V'r) f° r such a are displayed in Appendix ||. 

We can calculate the probability of observing the particle doing the quantum walk at any given point n = at 
from the wave function derived above. Below we give the asymptotic distribution for points a = n/t 
between —l/y/2 + e and l/y/2 — e, for an arbitrarily small constant e > 0. 

P(a,t) = |^ L (aM)| 2 + |^R(aM)| 2 
1 + (-l)^ 1 )* 

~ 1 — 77—, x 

7rt K 

(1 - af cos 2 {4>{a)t + vr/4) + (1 — a 2 ) cos 2 (4>(a) t + k a + tt/4) , (8) 
where 0(a) = — uJk a — ak a , and k a is the root of u>' k + a = in [0, n]. 

The (approximate) probability distribution P compares very well with simulation results even for small t, 
as is evident from Figure ||. The bias to the left in the probability distribution plotted in the figure is 
an artifact of the choice of initial chirality state (it was chosen to be "left"). If the particle begins in the 
chirality state ^(|L) + i|R)) , the distribution at any time turns out to be symmetric (cf. Section and 
Figure [5] in the appendix). Indeed, the Hadamard walk is an unbiased walk. 



3.2 Calculation of the moments 

We show in Appendix |C] that the net probability of the points n with a = n/t between —l/y/2 + e 
and l/y/2 — e, where e is an arbitrarily small constant, is 1 — — — — ^ , so the rest of the points do not 
contribute to any global properties of the distribution. Henceforth, we restrict ourselves to this interval. 

For the purposes of calculating the moments of the distribution, it will be convenient to decompose P as 

P{a,t) = P slow (a,t) + P fast (a,t), (9) 

where 

1 - a 

P s \ow{a,t) = „ (10) 

is a slowly varying (non-oscillating) function in a, and i-f as t is the remaining (quickly oscillating) component. 
By the analysis shown in Appendix |C|, we can show that any contribution to a moment from the "fast" 
component Pf ast is of lower order in t than the contribution from P s i w I n Figure |||, we compare P s low 
with P. 
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Figure 3: A comparison of the distributions P and P s \ ovl for t = 100. Only the probability at the even points 
is plotted, and P s iow is scaled by a factor of 2 because it has support on the odd points as well. 



The calculation of moments is further simplified by the following observation. Let p(a) = tP s \ ow (a,t). 
Then p is a probability density function over the interval between ^\j\[2. It is clearly non-negative, 
and we show below that it integrates to 1. This observation allows us to approximate the sums in the 
moment calculations by a Riemann integral, and the error so introduced is again a lower order term (cf. 
Appendix |C|) . 

To see that p(a) integrates to 1, note that ^(k a , a) = = — oj' ka — a, so 



Now, 




since the function cos k/y/l + cos 2 k is anti-symmetric around the point ir/2. 

The different moments for the density function p are now readily calculated by standard methods from 
complex residue theory. These are listed in Figure |4| for comparison with simulation results.^] This gives 
us the leading term for the moments for the distribution P. 



3.3 Mixing behaviour of the walk 

It is quite evident from Figure ^ that the probability distribution P is almost uniform over the interval 
between -1/^2, l/v^- We argue this formally below. 

Recall that the (5-mixing time (r^) of a random process is defined as the first time t such that the distribution 
at time t is at total variation distance (which is half the i\ distance) at most 5 from the uniform distribution. 

4 As mentioned before, the particle has a constant speed to the left, as indicated by its mean position, which is a result of 
its biased initial state. For an unbiased initial state, the mean would be zero. 
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0.293 



Figure 4: A table of moments calculated with the approximation by the density function p(a), which are 
compared with computer simulation results with t = 80. 



We claim that there is a constant 5 < 1 such that at time t, P is <5-close to the uniform distribution on the 
integer points between =pt/\/2- We emphasise that for classical symmetric random walk, the corresponding 
mixing time is quadratic in t. 

In order to show that the Hadamard walk is "mixed" at time t (in the sense described above), it suffices 
to show that a constant fraction [3 > of the points in the said interval have probability between c/y/2t 
and l/y/2t for a constant c > 0. A straightforward calculation shows that this implies that the i\ distance 
from uniform over 

[-t/y/2, t/y/2] is at most 2(1 - (5c), so that 5 = 1 - (3c < 1 is a constant. 

As in the previous section, we restrict ourselves to the interval where n/t = a £ [—l/y/2 + e, l/y/2 — e] 
(with e chosen to be a suitable constant) such that P(a,t) is at most \j\pit within the interval. Recall 
that the probability mass within this interval is 1 — — — which is a constant greater than (see 

Appendix ^| for a proof of this). Clearly, this cannot hold unless at least a constant fraction = 1 — — — c 

of the points within this interval have probability at least some constant < c < 1 over V2t. This 

completes the proof of the mixing nature of the Hadamard walk. 

4 Discussion 

In this paper, we defined and studied the behaviour of quantum walks on the line, and showed a linear 
time mixing behaviour for these walks. This has an immediate bearing on the similar walks defined on the 
circle. It is easy to derive the wave function for the walk on the circle, since it is exactly the wave function 
for the unbounded line wrapped around the circle. The linear time mixing for the circle is immediate. It 
is also possible to derive the wave function for an appropriately modified walk on the finite line using the 
solution to the unbounded walk. We expect that this walk be mixing in linear time as well, but leave the 
details to |l~6[| . (Note that because of boundary effects, the mixing of the walk on a finite line does not 
follow directly from our results if the walk begins close to the edges.) 

Aharonov et al. |]]] have defined quantum walks on general graphs in the same manner. It would be 
interesting to characterise the class of graphs on which quantum walks mix faster than classical ones. A 
promising candidate is the class of Cayley graphs, which seem to be amenable to analysis along the same 
lines as ours because of their rich group-theoretic structure. It would also be interesting to study how 
much speed-up is possible in general, and for different classes of graphs. The lower bound result of [U is a 
step in this direction. 

Our work leads to a plethora of questions regarding the properties of quantum walks. It is yet unclear which 
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of these are interesting from an algorithmic point of view. Any speed-up of a known classical algorithm 
based on random walk (such as those mentioned in Section [I]) seems to involve the analysis of quantum 
walk on graphs much more complex than the line. However, we believe such analysis is still tractable. 

Finally, we ask if it is possible to design quantum algorithms to generate certain desirable superpositions 
efficiently — those corresponding to distributions that are classically hard to sample from (much in the spirit 



of [15 1). It is conceivable that such algorithms follow an approach that is completely different from the 



standard Markov chain Monte Carlo method. 
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A Technical background 

In this section, we give details of the concepts and results from mathematical analysis that we use in the 
paper. 

A.l The Fourier transform 

Let /:2^Cbea complex valued function over the integers. Then its Fourier transform / : [— ir, ir] — ► C 
is defined as 

/>) = £/(n)e ifo \ 

n 

The corresponding inverse Fourier transform is then given by 

f(n) = ± £j~(k)e- ikn dk. 
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Note that this is the dual view of the Fourier transform on the space of functions on [— 7r,7r]. 

We will be concerned with functions / with finite support, i.e., functions which are zero except at finitely 
many points n. The Fourier transforms are extremely well behaved for this class of functions. For more 
details on the properties of these transforms, see, for example, 0). 



A. 2 Asymptotic expansion of integrals 

Studying the large time behaviour of quantum walks naturally leads us to consider the behaviour of integrals 
of the form 

J(f) = ( b g(k)e ltm dk (12) 

J a 

as t tends to infinity. There is a well-developed theory of the asymptotic expansion of integrals which 
allows us to determine, very precisely, the leading terms in the expansion of the integral in terms of simple 
functions of t (such as inverse powers of t). Below, we explain an important technique known as the method 
of stationary phase ||, |j] from this theory. 

Intuitively, the method may be understood as follows. The exponential in the integral is a rapidly oscillating 
function if t is large (and if (ft is n °t constant in any subinterval) . If g is a smooth function of k, then 
the contributions from adjacent subintervals nearly cancel each other out, and the major contribution to 
the value of the integral comes from the region where the oscillations are least rapid. The regions of slow 
oscillation occur precisely at the stationary points of the function (ft, i.e, points c where (ft'(c) = 0. (If no 
such point exists, we can get the asymptotic expansion for I in terms of inverse powers of t by repeated 
integration by parts, and the integral decays faster than j.) So the significant terms in the expansion of 
the integral come from a small interval around the stationary points. The "flatness" of (ft at a stationary 
point then determines the contribution of its neighbourhood to the integral. For example, if (ft"{c) ^ 0, 
then the integral goes as however if eft" {c) = but (ft"'(c) 0, then the integral goes as ^73. 

Without going into the details of its derivation, we state the leading term in the expansion of I(t) assuming 
that it has exactly one stationary point, and that it occurs at the left end point a of the interval. (Any 
integral may be written as a sum of such integrals.) We also assume that g is smooth and non-vanishing 
at a. Suppose a is a stationary point of order p—1, i.e., (ft' (a) = (ft" (a) = ■ ■ ■ = 0( p_1 )(a) = but (ft^ p \a) 7^ 0, 
then the dominant behaviour of I is given by 

I(t) ~ g(a) ^(a)±in/2p 

where we use the factor e i7r / 2 P (respectively, e~' m / 2 P) if <f>(p)( a ) > (if (ft^\a) < 0). For the reader interested 
in more details about this and other associated techniques, we recommend [||, ||. 



pi 



t\(ft(p){c 



r(i/p) 



p 



-00, 



(13) 



B Solution to the Hadamard walk 



The following are the two eigenvectors of the matrix for the Hadamard walk mentioned in Section 2J2 
(Recall that oJk is defined as the angle in [— |, ^] such that sin(wfc) = ^p) 



—j= ({I + cos 2 k) — cos k\J 1 + cos 2 k 



-ik 



V2i 



-ik 
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— j= ((I + cos 2 k) + cos 1 + cos 2 k 



-ik 



-ik 



-V2e luJk -e" 

We now get the expression for the Fourier transform at time t, given the initial state from the identity 



The expressions for ipL^R stated in Section [2.2| are derived by calculating this for a particle that begins 
the walk at the origin with chirality "left" (i.e., with ^f(k,0) = (1,0) T for all k). 

We now turn to determining the asymptotic behaviour of the wave function \& as t tends to +00. It suffices 
to consider an integral of the form 



I(a,t) 



7T ^7T 



(14) 



where g(k) is analytic, and a periodic function of period 2tt taken to be either even or odd, (f>(k, a) = 
— {oJk + oik), and a G [—1,1]. The integrals in the expressions for iph^R are exactly of this kind. 

It turns out that for large t, the integral I(a, t) has three distinct kinds of behaviour depending on the 
value of \a\, with sharp transitions from one behaviour to the other. It decays faster than any inverse 
polynomial in t for \a\ > + 0(i -2 / 3 ), goes as i -1 / 3 in an 0(i -2 / 3 ) interval around and as t -1 / 2 in 

the remaining interval. Below, we show these behaviours for / for \a\ > + e, \a\ = and |a| < — e, 
for any constant e > 0, and leave the more detailed analysis of the transitions from one behaviour to 



another to [16]. 

First, we concentrate on \a\ larger than by a constant. For this range of a, <f) does not have any 
stationary points, and we can use integration by parts to show the vanishing nature of I. 

Before we sketch a proof for this, we calculate some quantities that will also be of use later. 

cos k 



30 
dk 
d 2 cj> 
dk 2 
d 3 cj> 
W 



-(<4 + a) 
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y/l + cos 2 k 
sin A; 



a 



(1 + cos 2 A:) 3 / 2 
2 cos k(l + sin 2 k) 
(1 + cos 2 kf/ 2 



(15) 
(16) 
(17) 



Note that for \a\ > l/y/2 + e, ^ > e, so its inverse is an analytic function in k. Let 



u(k) 



11 dk 



and v(k) 



j<t>(k)t 



(In the above, we suppress the dependence of the functions on a for ease of notation.) Now, integrating 
by parts, 

TV \ f 71 dk dv r . . , f w dk du 

I(a,t) = j ^u- = [u{,)v{,)-u{-,)v(-,)]- j_^-v-. 

The first term above is because of periodicity, and the integrand in the second term is bounded (in 
magnitude) by c e /t for some constant c e that depends only on e. To prove a similar statement for any 
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greater power of t, we may use induction with the invariant that the integral resulting from the previous 
integration by parts is of the same form as I. 

We now look at the points a = l/\/2, — 1/\/2. At these points, (j> has a stationary point of order 2 
at k = it, 0, respectively, as may readily be verified. Using the method of stationary phase, we thus get the 
following leading term for / at these points: 



5Vrj 
3vr 



v / 2r(l/3) 



61 1 / 3 / 7T 7T 

- cos — =i H — 

-61 V3 
t 



Finally, we turn to the interval of most interest to us, [— l/y/2 + e, l/v2 — e]. When a lies in this region, 
4> has two stationary points k a ,—k a , where k a £ [0, w] and 



cos A; Q 



-a 



or 



The phase and its second derivative at the stationary points are: 

cp(±k a ,a) = T(u ka +ak a ) 

d 



^(±k a ,a) = ^Jl = ±{l-a 2 )Vl-2a 2 . 



(18) 

(19) 
(20) 



We abbreviate 4>(k a ,a) by 4>(a). We can again employ the method of stationary phase to get the dominant 
term in the expansion for /: 



I(a,t) 



27rfK' I 



2 cos((/>(a) t + vr/4) if g is even 
2isin(</>(a) t + vr/4) if g is odd 



(21) 



Using the above analysis, we can write out the asymptotic expressions for ifjh^R'- 

(1 - a) cos(cp(a) t + 7r/4) 



ip R (at,t) 



1 + (_!)(«+!)* 

> , X < 

2ntUl I 



a cos(<fi(a) t + vr/4) 
- Vl - 2a 2 sin(0(a) t + vr/4) 



(22) 



It is now straightforward to calculate the probability distribution P(a, t) induced on the line by observing 
the position of the particle after t steps of the walk. Equation |8| in Section [| gives the asymptotic expression 
for P(a,t). 

Our estimate of the probability distribution is compared with computer simulations in Figure ^ of Section [3|. 
In Figure [|, we give a plot of the distribution obtained with the symmetric initial chirality state 775 (|L) + 
i|R)). The figure illustrates the unbiased nature of the Hadamard walk. We leave the details of the 
calculation of the asymptotic behaviour of the distribution to |H| . 
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Figure 5: The probability distribution obtained from a computer simulation of the Hadamard walk with a 
symmetric initial condition. The number of steps in the walk was taken to be 100. Only the probability at 
the even points is plotted, since the odd points have probability zero. 



C Justification for approximations made in section |3| 

In calculating the moments of the probability distribution P(a, t) obtained by observing a particle doing 
the Hadamard walk, we made two approximations. First, we decomposed P into a slowly varying, non- 
oscillating part -P s iow an d a quickly oscillating part -Pf as t • We claimed that the contribution to the moments 
from Pf as t is of smaller order in t. Then, we approximated the moment sums involving P s i ow by a Reimann 
integral. In this section, we will sketch proofs that the error involved in the two approximations both go 
as 1/t times the quantity being estimated. 

As an important consequence of the error bound we derive, we show that the net probability of points in 
the range a G [—l/y/2 + e, l/y/2 — e] differs from the integral of the density function 



p(a) 



1 



a 



2c 



(23) 



over this interval by 0(l)/i. This integral equals 1 — — . 

In the following discussion, we assume that n/t,a £ [-l/y/2 + e, l/y/2-e] (th is range also applies to sums 
over n, integrals over a etc.). 

First, we bound the error in approximating a moment sum using P s i ow by an integral. We compute the 
moments in terms of fractions of the appropriate powers of t. For the mth moment (m > 0), 



£( w /*r p siow(n/ M ) 



a m p(a, t) da 



< E r +1)/t da \(n/t) m tP slow (n/t,t) - a m p(a)\ 

n Jn / t 

Kn+l)/t 

E / daA 1/t (a m p(a)) 

n Jn l t 



< 



(24) 



where 



A s f(a) = max |/(a + £)- /(a) | 
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It is easy to verify that the variation given by Aj/j for a m p(a) over a in the interval above is 0(1) /t. This 
gives an error bound of 0(l)/t in equation ( ^4|) above. 

We now show a bound of 0(1) /t on the error introduced by the first kind of approximation mentioned 
above, namely by ignoring Pf as t, the component of P with quickly oscillating factors. We bound one of the 
terms occurring in it; the others may be bounded similarly. Consider the following term obtained from 
equation (§) in Section ||, by expressing the sine and cosine functions in terms of exponentials: 

P(a,t) = a(1 ," 2 " ) e^ i , 

where 4>(a) = 4>(k a ,a). We estimate it's contribution to the mth moment. The basic observation is that 
each term in the moment sum such as a m P(a, t) is essentially the leading term in the asymptotic expansion 
of an integral of the type in equation (|il^) discussed Appendix [A|: 

-S 



a m P(a,t) ~ IT dkf(k)e 2i<t>{k ' a)t 
\Jt J6 



where 



2^(1 + 2^)^/4 



f(k) = — fc y fc v 

and 5 is a suitable constant. The expression for / above is derived using the identity a = —u>' k and 



comparing the term in the moment sum with the asymptotic form of the integral in equation (13) in 
Appendix [A]. The error introduced by this integral representation is 0(i~ 3 / 2 ). 

We may now bound the contribution to the moment sum due to P as 

-8 

2\u) k t—2\nk 



Y^(n/t) m P(a,t) ~ 4=Ef~ dk f( k )* 

Vt n J8 

J* 5 dkf(k)e~ 2iuJkt J2e 



1 

i r n-8 „-2ifc(/3t+l) _ p 2ifc/3t 

= Til dkm ^ h ' e-^-l ■ 

where = (l/y/2-e). We may now use the method of stationary phase again to derive the asymptotic 
form of the integral. As we expect from the analysis we did in Appendix [Bl the integral goes as 1/^Ji, and 
thus the entire expression goes as 1/t. 



D General walk on the line 

So far we have been describing the details of the Hadamard walk on the line. Here, we consider the 
properties of the general quantum walk on the line. As mentioned in Section |l|, it suffices to consider a 
particle that carries a "left" or "right" chirality, and at each time step undergoes a unitary transformation 
of its chirality state and then moves accordingly, one step to the right if it ends in "right" chirality state, 
and one step to the left if it ends in the "left" chirality state. The most general such unitary transformation 
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contains four parameters. However, it can be shown [|0| that if we focus on the probability distribution 
that results from these walks, then there is a single parameter that labels a family of walks, and every 
element of the family gives rise to identical probability distributions P{n, t) for the same initial conditions. 
The walks given by this parameter are also singled out for study in [p"3f. 



We therefore consider a representative element from each family, which we write as 

M(9) = To(I®U e ), where (25) 
T = T_ <g> |L)(L| + T+ <g> |R)(R|, and (26) 
U e = e 1 ^. (27) 

The operator M(9) evolves the state by one time step, i.e., + 1)) = M(9)\^f(t)}, the operator T 
implements the left/right shift {T + \n) = |n+l) while T_|n) = \n— 1)), and Ug is the unitary transformation 
acting on the chirality state. The Hadamard walk belongs to the family of walks parametrised by 9 = \. 
To show this we write the Hadamard transform as 

H = ~^(\ _\ ^ = -ie^^e'i^ = -i£? a 'U%. (28) 

Thus, by a suitable redefinition of phase of the chirality states at time t + 1, the additional rotation that 
multiplies Uil can be absorbed. Clearly, the probability distributions for finding a particle at a point are 
unaffected by this redefinition of phase. Therefore, the Hadamard walk is equivalent to the walk with 9 = ^. 

The behaviour of the walk with 9 = and 9 = n are easily analysed. When 9 = 0, the unitary transforma- 
tion rotating the chirality states is the identity, and hence a particle starting with "left" ( "right" ) chirality, 
after t time steps, moves t units to the left (right). This walk is clearly not mixing, since the probability 
is concentrated at one or two points at any given instant, and resembles more the biased classical walk 
(though we wish to emphasise here that all the quantum walks we consider are unbiased, a point which we 
return to below). Now, when 9 = ir, the unitary transformation flips the chirality at each time step. Hence 
a particle zig zags between the origin and a neighbouring point, and thus has a non-zero probability on a 
maximum of three points. This walk too is clearly not mixing. 

For intermediate values of 9 the behaviour is more involved, as our analysis of the Hadamard walk indicates. 
In fact, the qualitative behaviour of the Hadamard walk is generic to all 9 excluding the singular values 0, 7r 
that we discussed above. All of these walks are mixing, and mix in linear time. 

We can calculate very easily some properties of the general walk, following the analysis for the Hadamard 
walk. We first note that the eigenvalues of Mk{9) for the general walk given by 9 are 

\\{9),\l{9) = where (29) 

Q 

cos{oJk(9)) = cos -cos A: (30) 

where u>k(6) £ [0, 7r]. 

It is possible to deduce the relevant component P s i ow of the asymptotic probability distribution P from 
the eigenvalues in equation (|2^), as is evident from the case of the Hadamard walk. The distribution P s i ow 
beginning with symmetric initial conditions is: 

W<M) = ^r^y (31) 
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and has support on the interval a S (— cos|, cos|). Here k a is given by the usual stationary point 
condition: 

doj k (9) 



dk 

and the derivatives of 0J k (9) with respect to k are 



= a, (32) 



, , , sin k cos o , , 

<4 = ■ , * 33 
sm.(uk{0)) 

J k \9) = cot(u k (9))[l-(uj' k (e)n (34) 



Thus the stationary points are k a ,ir — k a , where k a £ [— n/2, tt/2] and 



• 7 & /or-\ 

smk a = tan - . 35 



The second derivative of the phase at the stationary points is (9) where 



"1(9) = i-i/cos^- a 2. (36) 



sin - 2 



We can readily verify that the probability sums to one, i.e., 

r f 1 f 77 / 2 dk 

t daP Aow (a,t) = da = / — - = 1, (37) 

J J TXUl a {9) J -ir/2 TT 

where we have used the fact that u/ ka = Thus, as in the case of the Hadamard walk, almost all 

the probability is confined to the interval between =p cos | and is close to uniform over it. However, the 
"height" of the distribution is proportional to sin |, and the distance from uniform increases as 9 — > 0. On 
the other hand, as 9 — > it, the "width" of the distribution (which is proportional to 2 cos |) goes to zero, 
again increasing the distance from uniform. 

The distribution P s i ow can, as before, be used to calculate moments. For example, the mean position of 
the particle (\a\) is given by 

dk 9 
da\a\tP slow (a,t) = 2 — J ka {9) = 1--, (38) 

JO TT TT 

which matches the expected value at the points 9 = 0, 5,7r. 

Finally we comment on the unbiased nature of the quantum walks defined in this manner. An unbiased 
walk (starting at the origin) does not distinguish between the left and right halves of the line, i.e., it is 
reflection symmetric. In order to capture this into a condition on unbiased walks, it is necessary to also 
account for the initial chirality state of the particle which can lead to a biased probability distribution, 
although the rules governing the walk are intrinsically reflection symmetric. Therefore, although such a 
symmetric walk can produce a biased probability distribution with some particular initial condition, there 
exists another initial condition for which the bias is exactly reversed. In more formal terms, if we describe 
the quantum walk by the matrices M k in equation (|l|) in Section ^, then we need to find a unitary matrix 
S such that: 

&M k S = M_ k . (39) 
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The unitary matrix S performs a rotation on the chirality state that reverses any bias arising from the 
choice of initial condition. If no such matrix exists, then the walk intrinsically distinguishes left from right 
and is hence not an unbiased walkQ. An immediate consequence is that if we start with an initial chirality 
state that is an eigenvector of S, then a symmetric probability distribution is generated. For the family 
of walks just considered, the unitary matrix S = a y will suffice since it only interchanges T + and T_ 
and does nothing else. Therefore, starting in the initial chirality state ^(0, ±i) T generates a symmetric 
distribution, as shown for the Hadamard walk in Figure [|. 



5 Actually the condition is a little more general, an overall sign change can also be permitted, i.e., S' M^S = ±M_fc. In 
fact, the minus sign appears for the Hadamard walk, where S = a y . 
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